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Abstract. These lectures notes aim at introducing Levy processes in 
an informal and intuitive way, accessible to non-specialists in the field. 
In the first part, we focus on the theory of Levy processes. We analyze 
p^ a 'toy' example of a Levy process, viz. a Levy jump-diffusion, which yet 

^vi offers significant insight into the distributional and path structure of a 

Levy process. Then, we present several important results about Levy 
processes, such as infinite divisibility and the Levy-Khintchine formula, 
the Levy-Ito decomposition, the Ito formula for Levy processes and Gir- 
^H sanov's transformation. Some (sketches of) proofs are presented, still 

r-f\ the majority of proofs is omitted and the reader is referred to textbooks 

instead. In the second part, we turn our attention to the applications 

of Levy processes in financial modeling and option pricing. We discuss 

Ph how the price process of an asset can be modeled using Levy processes 

Ph and give a brief account of market incompleteness. Popular models in 

(-H the literature are presented and revisited from the point of view of Levy 

■ H processes, and we also discuss three methods for pricing financial deriva- 

I fives. Finally, some indicative evidence from applications to market data 

O"* is presented. 
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Part 1. Theory 



1. Introduction 



Levy processes play a central role in several fields of science, such as 
physics, in the study of turbulence, laser cooling and in quantum field theory; 
in engineering, for the study of networks, queues and dams; in economics, for 
continuous time-series models; in the actuarial science, for the calculation 
of insurance and re-insurance risk; and, of course, in mathematical finance. 
A comprehensive overview of several applications of Levy processes can be 
found in Prabhu (1998), in Barndorff-Nielsen, Mikosch, and Resnick (2001), 
in Kyprianou, Schoutens, and Wilmott (2005) and in Kyprianou (2006). 
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Figure 1.1. USD/JPY exchange rate, Oct. 1997-Oct. 2004. 



In mathematical finance, Levy processes are becoming extremely fashion- 
able because they can describe the observed reality of financial markets in 
a more accurate way than models based on Brownian motion. In the 'real' 
world, we observe that asset price processes have jumps or spikes, and risk 
managers have to take them into consideration; in Figure 1.1 we can observe 
some big price changes (jumps) even on the very liquid USD/JPY exchange 
rate. Moreover, the empirical distribution of asset returns exhibits fat tails 
and skewness, behavior that deviates from normality; see Figure 1.2 for a 
characteristic picture. Hence, models that accurately fit return distributions 
are essential for the estimation of profit and loss (P&L) distributions. Simi- 
larly, in the 'risk-neutral' world, we observe that implied volatilities are con- 
stant neither across strike nor across maturities as stipulated by the Black 
and Scholes (1973) (actually, Samuelson 1965) model; Figure 1.3 depicts a 
typical volatility surface. Therefore, traders need models that can capture 
the behavior of the implied volatility smiles more accurately, in order to 
handle the risk of trades. Levy processes provide us with the appropriate 
tools to adequately and consistently describe all these observations, both in 
the 'real' and in the 'risk-neutral' world. 

The main aim of these lecture notes is to provide an accessible overview 
of the field of Levy processes and their applications in mathematical finance 
to the non-specialist reader. To serve that purpose, we have avoided most 
of the proofs and only sketch a number of proofs, especially when they offer 
some important insight to the reader. Moreover, we have put emphasis on 
the intuitive understanding of the material, through several pictures and 
simulations. 

We begin with the definition of a Levy process and some known exam- 
ples. Using these as the reference point, we construct and study a Levy 
jump-diffusion; despite its simple nature, it offers significant insights and an 
intuitive understanding of general Levy processes. We then discuss infinitely 
divisible distributions and present the celebrated Levy-Khintchine formula, 
which links processes to distributions. The opposite way, from distributions 
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Figure 1.2. Empirical distribution of daily log-returns for 
the GBP/USD exchange rate and fitted Normal distribution. 



to processes, is the subject of the Levy-Ito decomposition of a Levy pro- 
cess. The Levy measure, which is responsible for the richness of the class of 
Levy processes, is studied in some detail and we use it to draw some conclu- 
sions about the path and moment properties of a Levy process. In the next 
section, we look into several subclasses that have attracted special atten- 
tion and then present some important results from semimartingale theory. 
A study of martingale properties of Levy processes and the Ito formula for 
Levy processes follows. The change of probability measure and Girsanov's 
theorem are studied is some detail and we also give a complete proof in the 
case of the Esscher transform. Next, we outline three ways for constructing 
new Levy processes and the first part closes with an account on simulation 
methods for some Levy processes. 

The second part of the notes is devoted to the applications of Levy pro- 
cesses in mathematical finance. We describe the possible approaches in mod- 
eling the price process of a financial asset using Levy processes under the 
'real' and the 'risk-neutral' world, and give a brief account of market incom- 
pleteness which links the two worlds. Then, we present a primer of popular 
Levy models in the mathematical finance literature, listing some of their 
key properties, such as the characteristic function, moments and densities 
(if known). In the next section, we give an overview of three methods for pric- 
ing options in Levy-driven models, viz. transform, partial integro-differential 
equation (PIDE) and Monte Carlo methods. Finally, we present some em- 
pirical results from the application of Levy processes to real market financial 
data. The appendices collect some results about Poisson random variables 
and processes, explain some notation and provide information and links re- 
garding the data sets used. 

Naturally, there is a number of sources that the interested reader should 
consult in order to deepen his knowledge and understanding of Levy pro- 
cesses. We mention here the books of Bertoin (1996), Sato (1999), Apple- 
baum (2004), Kyprianou (2006) on various aspects of Levy processes. Cont 
and Tankov (2003) and Schoutens (2003) focus on the applications of Levy 
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Figure 1.3. Implied volatilities of vanilla options on the 
EUR/USD exchange rate on November 5, 2001. 



processes in finance. The books of Jacod and Shiryaev (2003) and Prot- 
ter (2004) are essential readings for semimartingale theory, while Shiryaev 
(1999) blends semimartingale theory and applications to finance in an im- 
pressive manner. Other interesting and inspiring sources are the papers by 
Eberlein (2001), Cont (2001), Barndorff-Nielsen and Prause (2001), Carr et 
al. (2002), Eberlein and 6zkan(2003) and Eberlein (2007). 



2. Definition 

Let (Jl,^", F,P) be a filtered probability space, where T = Tt and the 
filtration F = (^ r t)ie[o,r] satisfies the usual conditions. Let T G [0,oo] denote 
the time horizon which, in general, can be infinite. 

Definition 2.1. A cadlag, adapted, real valued stochastic process L = 
(L t )o<t<T with Lq = a.s. is called a Levy process if the following con- 
ditions are satisfied: 

(LI): L has independent increments, i.e. Lf — L s is independent of J- s 
for any < s < t < T. 

(L2): L has stationary increments, i.e. for any < s,t < T the distri- 
bution of L t + S — L t does not depend on t. 

(L3): L is stochastically continuous, i.e. for every < t < T and e > 0: 



linu 



P L t -LJ >e 



0. 



The simplest Levy process is the linear drift, a deterministic process. 
Brownian motion is the only (non-deterministic) Levy process with contin- 
uous sample paths. Other examples of Levy processes are the Poisson and 
compound Poisson processes. Notice that the sum of a linear drift, a Brow- 
nian motion and a compound Poisson process is again a Levy process; it 
is often called a "jump-diffusion" process. We shall call it a u Levy jump- 
diffusion" process, since there exist jump-diffusion processes which are not 
Levy processes. 
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Figure 2.4. Examples of Levy processes: linear drift (left) 
and Brownian motion. 




Figure 2.5. Examples of Levy processes: compound Poisson 
process (left) and Levy jump-diffusion. 



3. 'Toy' example: a Levy jump-diffusion 

Assume that the process L = (Lt)o<t<T is a Levy jump- diffusion, i.e. a 
Brownian motion plus a compensated compound Poisson process. The paths 
of this process can be described by 



(3.1) 



Nt 
L t = bt + aW t + ( X) Jfe _ tXK ) 



fc=i 



where b G R, a £ M^o 5 W = (Wt)o<t<T is a standard Brownian motion, 
N = (N t )o<t<T is a Poisson process with parameter A (i.e. E[iVj] = At) 
and J = (Jfe)fe>i is an i.i.d. sequence of random variables with probability 
distribution F and E[J] = re < cxd. Hence, F describes the distribution of 
the jumps, which arrive according to the Poisson process. All sources of 
randomness are mutually independent. 

It is well known that Brownian motion is a martingale; moreover, the 
compensated compound Poisson process is a martingale. Therefore, L = 
(Lt)o<t<T is a martingale if and only if b = 0. 
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The characteristic function of Lt is 
JE[e iuLt ] = E exp (iu(bt + aW t + ^ J k - tXn)) 



fc=i 



N t 



exp [i-u&iJE exp (iuaWt) exp (iu( \^ </& — £Ak)) 



fc=i 



since all the sources of randomness are independent, we get 



exp [i-ufeiJE 



exp (iuaWt 



E 



N t 



exp (iu y^ Jk — iutXt 



fc=i 



taking into account that 

E[ e iwrW *] = e"^ 2 " 2 *, W t ~ Normal(0, t) 

E[e iu ^=i Jfc ] = e^^I^- 7 - 1 ]), iV t ~ Poisson(At) 
(cf. also Appendix B) we get 



exp [iubt] exp 
exp [iubt] exp 



2 2. 
-u a t 



2 2. 

u a t 



exp 
exp 



At(E[e iuJ - 1] - wE[J]) 
At(E[e iuJ -l-inJ]) 



and because the distribution of J is -F we have 

1 



exp [iubt] exp 



2 2. 

-u a t 



exp 



At / (e iu *-l-iux)F(dx) 



Now, since t is a common factor, we re-write the above equation as 



(3.2) lE[e iuLt ] =exp 



t ( iub 



..v 



+/<*-- 1-*,)^: 



Since the characteristic function of a random variable determines its dis- 
tribution, we have a "characterization" of the distribution of the random 
variables underlying the Levy jump-diffusion. We will soon see that this dis- 
tribution belongs to the class of infinitely divisible distributions and that 
equation (3.2) is a special case of the celebrated Levy-Khintchine formula. 

Remark 3.1. Note that time factorizes out, and the drift, diffusion and 
jumps parts are separated; moreover, the jump part factorizes to expected 
number of jumps (A) and distribution of jump size (F). It is only natural 
to ask if these features are preserved for all Levy processes. The answer is 
yes for the first two questions, but jumps cannot be always separated into a 
product of the form A x F. 
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4. Infinitely divisible distributions and the Levy-Khintchine 

FORMULA 

There is a strong interplay between Levy processes and infinitely divisible 
distributions. We first define infinitely divisible distributions and give some 
examples, and then describe their relationship to Levy processes. 

Let X be a real valued random variable, denote its characteristic function 
by ipx and its law by Px, hence <px(u) = j^ elux Px(dx). Let \x * v denote 
the convolution of the measures \x and v, i.e. (/x * v)(A) = j K v{A — x)/J,(dx). 

Definition 4.1. The law Px of a random variable X is infinitely divisible, 
if for all n G N there exist i.i.d. random variables X[ , . . . ,Xn such 
that 

(4.1) 



x d (i/„) + _ + x ( 1/n) _ 



Equivalently, the law Px of a random variable X is infinitely divisible if for 
all n G N there exists another law P X (i/„) of a random variable X^ ' n > such 
that 



(4.2) 



Px = P 



xC 1 /™) 



* ... * P 



XO-/"-) 



Alternatively, we can characterize an infinitely divisible random variable 
using its characteristic function. 

Characterization 4.2. The law of a random variable X is infinitely divis- 
ible, if for all n G N, there exists a random variable X^ l ' n > , such that 



(4.3) 



vx(u) = yp x ^/^( u )) 



Example 4.3 (Normal distribution). Using the characterization above, we 
can easily deduce that the Normal distribution is infinitely divisible. Let 
X ~ Normal^, a 2 ), then we have 

1 



ipx(u) = exp 



lUfl 



exp 



iu- 



2 

t - 

n 



uV" 


= exp 


1 2 
— u 

2 


a 2 i 
n . 


r 



n[iu 



fi 



n 



1 2<A 
-U — 

2 n ' 



<p X (i/n)(u) 



where X^ 1 ' n > ~ Normal( 



/J; CT^ 

n ' n 



Example 4.4 (Poisson distribution). Following the same procedure, we can 
easily conclude that the Poisson distribution is also infinitely divisible. Let 
X ~ Poisson(A), then we have 



<fx(u) = exp A(e m - 1 



rA 



exp 



n 



(e* 



i: 



^ x( i/n)(n) 
where A^ 1 /™) ~ Poisson(^). 
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Remark 4.5. Other examples of infinitely divisible distributions are the 
compound Poisson distribution, the exponential, the T-distribution, the geo- 
metric, the negative binomial, the Cauchy distribution and the strictly stable 
distribution. Counter-examples are the uniform and binomial distributions. 

The next theorem provides a complete characterization of random vari- 
ables with infinitely divisible distributions via their characteristic functions; 
this is the celebrated Levy-Khintchine formula. We will use the following 
preparatory result (cf. Sato 1999, Lemma 7.8). 

Lemma 4.6. If (Pk)k>o is a sequence of infinitely divisible laws and P^ — ► 
P, then P is also infinitely divisible. 

Theorem 4.7. The law Px of a random variable X is infinitely divisible if 
and only if there exists a triplet (6, c, u), with b £ M, c G M^o & n d a measure 
satisfying ^({0}) = and j R (l A |x| 2 )z/(dx) < oo, such that 

(4.4) TE[e iuX ] = exp \ibu - V + / ^ ~ l ~ iua:1 {M<i}M da; ) 

K 

Sketch of Proof. Here we describe the proof of the "if part, for the full proof 
see Theorem 8.1 in Sato (1999). Let (£ n )neN be a sequence in M, monotonic 
and decreasing to zero. Define for all u £ M and n £ N 

fx n {u)=exp iu(b- / xv(dx)\ - — + / (e lux - l)u(dx) . 

£ n <|a;|<l |a;|>en 

Each (px n is the convolution of a normal and a compound Poisson distri- 
bution, hence cpx n is the characteristic function of an infinitely divisible 
probability measure Px n ■ We clearly have that 

lim ip Xn (u) =ip x (u); 

then, by Levy's continuity theorem and Lemma 4.6, (px is the characteristic 
function of an infinitely divisible law, provided that ipx is continuous at 0. 
Now, continuity of (px at boils down to the continuity of the integral 
term, i.e. 



i/> v (u) = / (e iux - 1 - iuxl { \ xl<1} )u(dx) 



(e lux - 1 - iux)u{dx) + / {e tux - l)i/(dx). 

{|x|<l} {|x|>l} 

Using Taylor's expansion, the Cauchy-Schwarz inequality, the definition of 
the Levy measure and dominated convergence, we get 

\Mu)\<l I \u 2 x 2 \u(dx)+ J \e iux - l|i/(dx) 

{\x\<l} {\x\>l} 

I |2 /• p 

<^ / |x 2 |i/(dar)+ / \e iux - l|i/(dx) 

{\x\<l} {\x\>l} 

— >0 as u^O. □ 
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The triplet (b, c, v) is called the Levy or characteristic triplet and the 
exponent in (4.4) 

xi c f 
(4.5) ip(u) =iub- — + / {e lux - 1 - iuxl { \ x \ <1} )v(6x) 



is called the Levy or characteristic exponent. Moreover, b £ R is called the 
drift term, c G M^o the Gaussian or diffusion coefficient and i/ the Levy 
measure. 

Remark 4.8. Comparing equations (3.2) and (4.4), we can immediately 
deduce that the random variable Lt of the Levy jump- diffusion is infinitely 
divisible with Levy triplet b = b • t,c = a 2 • t and v = (XF) • t. 

Now, consider a Levy process L = {Lt)o<t<T] for any n € N and any 
< t < T we trivially have that 

(4.6) L t = Lt+ {Lit - Lt) + . . . + {L t - L (n - 1)t ). 

n n n n 

The stationarity and independence of the increments yield that {Ltk — 

n 

L t (k-i))k>i is an i.i.d. sequence of random variables, hence we can conclude 

n 

that the random variable Lt is infinitely divisible. 

Theorem 4.9. For every Levy process L = {Lt)o<t<T, we have that 

(4.7) E[e iuLi ] = e^ (u) 



exp 



u 2 c 



2 



t(ibu --=-+/ {e tux - 1 - iwxl { | x |<i})i/(dx)) 



where ijj{u) is the characteristic exponent of L\, a random variable with an 
infinitely divisible distribution. 

Sketch of Proof. Define the function cf) u {t) = (pL t {u), then we have 

(4.8) (t> u {t + s)= TE[e iuLt+s ] = ^<L t+s -L s ) e iuL^ 

= E[e iu(L '+^ Ls) ]E[e iuLs ] = 0„(i)#u(*)- 

Now, (j) u {0) = 1 and the map t h^ <p u {t) is continuous (by stochastic con- 
tinuity). However, the unique continuous solution of the Cauchy functional 
equation (4.8) is 

(4.9) M (t)= e M ("), where : R -► C. 

Since L\ is an infinitely divisible random variable, the statement follows. □ 

We have seen so far, that every Levy process can be associated with the 
law of an infinitely divisible distribution. The opposite, i.e. that given any 
random variable X, whose law is infinitely divisible, we can construct a Levy 
process L = {L t )o<t<T such that C{L\) := £{X), is also true. This will be 
the subject of the Levy-Ito decomposition. We prepare this result with an 
analysis of the jumps of a Levy process and the introduction of Poisson 
random measures. 
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5. Analysis of jumps and Poisson random measures 

The jump process AL = (ALt)o<t<T associated to the Levy process L is 
defined, for each < t < T, via 

AL t = L t -L t _, 

where Lt- = lim s ^ t L s . The condition of stochastic continuity of a Levy 
process yields immediately that for any Levy process L and any fixed t > 0, 
then ALt = a.s.; hence, a Levy process has no fixed times of discontinuity. 
In general, the sum of the jumps of a Levy process does not converge, in 
other words it is possible that 



£ial s 



but we always have that 



oo a.s. 

s<t 



\] | Ai s | 2 < oo a.s. 



s<t 

which allows us to handle Levy processes by martingale techniques. 

A convenient tool for analyzing the jumps of a Levy process is the random 
measure of jumps of the process. Consider a set A 6 23(M\{0}) such that 
^ A and let < t < T; define the random measure of the jumps of the 
process L by 

(5.1) /z L (w; t, A) = #{0 < s < t; AL s (uj) £ A} 

= J2 U(AL»); 

s<t 

hence, the measure /j, (u>; t, A) counts the jumps of the process L of size in 
A up to time t. Now, we can check that n L has the following properties: 

fi L (t, A) - fi L (s, A) € o-{{L u - L v ;s <v < u < t}) 

hence \i (t, A) — n (s, A) is independent of J- s , i.e. (/,(-, A) has independent 
increments. Moreover, \x (t,A) — fi (s,A) equals the number of jumps of 
L s + U — L s in A for < u < t — s; hence, by the stationarity of the increments 
of L, we conclude: 

£(n L (t, A) - n L {s, A)) = £(/i L (t - s, A)) 

i.e. fj, L (-,A) has stationary increments. 

Hence, // (■, A) is a Poisson process and \x is a Poisson random measure. 
The intensity of this Poisson process is v(A) = E[/i (1, A)]. 

Theorem 5.1. The set function A h^ h l {uj; t, A) defines a a-finite measure 
on M\{0} for each (u>,t). The set function v(A) = E[^ L (1, J 4)] defines a 
a-finite measure on M\{0}. 

Proof. The set function A h^ /j, L (uj;t, A) is simply a counting measure on 
6(M\{0}); hence, 

TE[fi L (t,A)}= I fi L (uj;t,A)dP(uj) 

is a Borel measure on i3(M\{0}). □ 
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Definition 5.2. The measure v denned by 

v{A) = E[^ L (1,^)] = JE[J2 MAL»)] 



S<1 



is the Levy measure of the Levy process L. 

Now, using that [i (t,A) is a counting measure we can define an in- 
tegral with respect to the Poisson random measure \x L '. Consider a set 
A G #(M\{0}) such that ^ A and a function / : R -> R, Borel measur- 
able and finite on A. Then, the integral with respect to a Poisson random 
measure is defined as follows: 

(5.2) J f(x)fi L (Lo; t, dx) = J2 f(AL s )l A (AL s (u)). 

A s ^ 

Note that each f A f(x)fi L (t,dx) is a real-valued random variable and gen- 
erates a cadlag stochastic process. We will denote the stochastic process by 

Jo Ja /( x )^ L ( ds > dx ) = (Jo Ia /( x )^ L ( ds ' dx)) <t<T- 

Theorem 5.3. Consider a set A £ B(M\{0}) with ^ A and a function 
f : R — > W, Borel measurable and finite on A. 

A. The process (L f A f(x)/i L (ds,dx))o<t<T is a compound Poisson process 
with characteristic function 



(5.3) E expfiw / / f(x)^(ds, dx) 

B. If f € L l (A), then 
ft 

(5.4) E 
/o 

C. Iff€ L 2 (A), then 
ft 

(5.5) Var ' 



exp ( t I (e 



iuf(x) 



l)u(dx] 



f(x)^(ds,dx) =t / f(x)u(dx). 



f(x)ii L (ds,dx) 



t I \f(x)\ 2 v(dx). 
IA 



Sketch of Proof. The structure of the proof is to start with simple functions 
and pass to positive measurable functions, then take limits and use domi- 
nated convergence; cf. Theorem 2.3.8 in Applebaum (2004). □ 

6. The Levy-Ito decomposition 

Theorem 6.1. Consider a triplet (b,c,is) where b G R, c £ R^o an d v is a 
measure satisfying ^({0}) = and j K (l A\x\ 2 )v(dx) < oo. Then, there exists 
a probability space (O, J 7 , P) on which four independent Levy processes exist, 
ZA 1 ', u- 2 >, L^> and L^>, where ZA 1 ) is a constant drift, L^ 2 > is a Brownian 
motion, L^> is a compound Poisson process and L^> is a square integrable 
(pure jump) martingale with an a.s. countable number of jumps of magnitude 
less than 1 on each finite time interval. Taking L = Z/ 1 ' + L^ 2 > + L^> + L^> , 
we have that there exists a probability space on which a Levy process L = 
(Lt)o<t<T with characteristic exponent 



(6.1) 



ip(u) = iub - 



u 2 c 



+ /(e 



1 — iuxl 



{\x\<l 



} )i/(da 
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for all tiGl, is defined. 

Proof. See chapter 4 in Sato (1999) or chapter 2 in Kyprianou (2006). □ 

The Levy-Ito decomposition is a hard mathematical result to prove; here, 
we go through some steps of the proof because it reveals much about the 
structure of the paths of a Levy process. We split the Levy exponent (6.1) 
into four parts 

y, = y,« + y/ 2 ) + ^) + y,( 4 ) 

where 

i; (3 \u)= I (e iux -l)u(dx), 

\x\>l 

^\u) = J (e iux - 1 - iux)v{dx). 

\x\<l 

The first part corresponds to a deterministic linear process (drift) with pa- 
rameter b, the second one to a Brownian motion with coefficient \fc and 
the third part corresponds to a compound Poisson process with arrival rate 
A := i/(R\(-l, 1)) and jump magnitude F(dx) := ^^j) 1 {\x\>i}- 

The last part is the most difficult to handle; let AZ/ 4 ) denote the jumps 
of the Levy process Lft\ that is AL t = L t — L t J, and let /x L denote 
the random measure counting the jumps of L' 4 \ Next, one constructs a 
compensated compound Poisson process 

L ^= E A 4% >|AL ( 4 » l>£} -< / ww 

0<S<t -wl 1^ 

- - l>|x|>e 

t 

x^\dxAs)-t{ J xu(dx)) 

l>|x|>e l>|x|>e 

and shows that the jumps of L^> form a Poisson process; using Theorem 
5.3 we get that the characteristic exponent of IA 4 ' e ) is 

**><„)= / («»-.-)*). 

e<|x|<l 

Then, there exists a Levy process L^> which is a square integrable martingale, 
such that Z/ 4 ' e ) — ► L^ 4 ' uniformly on [0, T] as e — > 0+. Clearly, the Levy 
exponent of the latter Levy process is ip^> . 

Therefore, we can decompose any Levy process into four independent 
Levy processes L = L^> + L^ + L^ + L^\ as follows 

t t 

(6.2) L t = bt + y/cW t + J J xfi L (ds,dx)+ J f x(fi L - u L )(ds,dx) 

| X |>1 |x|<l 
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Figure 7.6. The distribution function of the Levy measure 
of the standard Poisson process (left) and the density of the 
Levy measure of a compound Poisson process with double- 
exponentially distributed jumps. 



Figure 7.7. The density of the Levy measure of an NIG 
(left) and an a-stable process. 



where v (ds,dx) = u(dx)ds. Here L"> is a constant drift, L' 2 ) a Brownian 
motion, L^> a compound Poisson process and Is-' a pure jump martingale. 
This result is the celebrated Levy-Ito decomposition of a Levy process. 

7. The Levy measure, path and moment properties 
The Levy measure v is a measure on R that satisfies 

(7.1) i/({0}) = and / (1 A \x\ 2 )v{dx) < oo. 



Intuitively speaking, the Levy measure describes the expected number of 
jumps of a certain height in a time interval of length 1. The Levy measure 
has no mass at the origin, while singularities (i.e. infinitely many jumps) can 
occur around the origin (i.e. small jumps). Moreover, the mass away from 
the origin is bounded (i.e. only a finite number of big jumps can occur). 

Recall the example of the Levy jump-diffusion; the Levy measure is v(dx) = 
A x F(dx); from that we can deduce that the expected number of jumps is 
A and the jump size is distributed according to F. 

More generally, if v is a finite measure, i.e. A := v(M) = J" R v{dx) < oo, 

then we can define F(dx) := \ , which is a probability measure. Thus, 
A is the expected number of jumps and F(dx) the distribution of the jump 
size x. If v (R) = oo, then an infinite number of (small) jumps is expected. 
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Figure 7.8. The Levy measure must integrate \x\ 2 A 1 (red 
line); it has finite variation if it integrates \x\ A 1 (blue line); 
it is finite if it integrates 1 (orange line) . 

The Levy measure is responsible for the richness of the class of Levy 
processes and carries useful information about the structure of the process. 
Path properties can be read from the Levy measure: for example, Figures 
7.6 and 7.7 reveal that the compound Poisson process has a finite number 
of jumps on every time interval, while the NIG and a-stable processes have 
an infinite one; we then speak of an infinite activity Levy process. 

Proposition 7.1. Let L be a Levy process with triplet (b,c,v). 

(1) Lf v(M.) < oo ; then almost all paths of L have a finite number of 
jumps on every compact interval. Ln that case, the Levy process has 
finite activity. 

(2) If u(R) = oo ; then almost all paths of L have an infinite number of 
jumps on every compact interval. In that case, the Levy process has 
infinite activity. 

Proof. See Theorem 21.3 in Sato (1999). □ 

Whether a Levy process has finite variation or not also depends on the 
Levy measure (and on the presence or absence of a Brownian part). 

Proposition 7.2. Let L be a Levy process with triplet (b,c,v). 

(1) If c = and Ji x i < i \x\v(dx) < oo, then almost all paths of L have 
finite variation. 

(2) If c t^ or J\ x \ <1 \x\v(dx) = oo, then almost all paths of L have 
infinite variation. 

Proof. See Theorem 21.9 in Sato (1999). □ 

The different functions a Levy measure has to integrate in order to have 
finite activity or variation, are graphically exhibited in Figure 7.8. The com- 
pound Poisson process has finite measure, hence it has finite variation as 
well; on the contrary, the NIG Levy process has an infinite measure and has 
infinite variation. In addition, the CGMY Levy process for < Y < 1 has 
infinite activity, but the paths have finite variation. 
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The Levy measure also carries information about the finiteness of the 
moments of a Levy process. This is particularly useful information in math- 
ematical finance, related to the existence of a martingale measure. 

The finiteness of the moments of a Levy process is related to the finiteness 
of an integral over the Levy measure (more precisely, the restriction of the 
Levy measure to jumps larger than 1 in absolute value, i.e. big jumps). 

Proposition 7.3. Let L be a Levy process with triplet (b,c,v). Then 

(1) Lt has finite p-th moment for p £ M^o (E|L(| P < oo) if and only if 

I\x\>i MMdz) < oo. 

(2) Lt has finite p-th exponential moment for p£R (E[e p *] < oo) if 
and only if j\, >l e px v{dx) < oo. 

Proof. The proof of these results can be found in Theorem 25.3 in Sato 
(1999). Actually, the conclusion of this theorem holds for the general class of 
submultiplicative functions (cf. Definition 25.1 in Sato 1999), which contains 
exp(px) and \x\ p V 1 as special cases. □ 

In order to gain some understanding of this result and because it blends 
beautifully with the Levy-Ito decomposition, we will give a rough proof of 
the sufficiency for the second statement (inspired by Kyprianou 2006). 

Recall from the Levy-Ito decomposition, that the characteristic exponent 
of a Levy process was split into four independent parts, the third of which 
is a compound Poisson process with arrival rate A := i/(R\(— 1,1)) and 
jump magnitude F(dx) := aruJ'i iii -*-{M>i}- Finiteness of E[e pLt ] implies 



finiteness of E|e pL « 1, where 



E^Ve^V^f- I / e px F(dx) 




-^ST— I / „/'■'■ 



^ k\ 

k>0 



c 



l{|x|>ip(dz) 



Since all the summands must be finite, the one corresponding to k = 1 must 
also be finite, therefore 

xt J e pz 



: l{| :E |> 1 }^(dx) < oo =>■ / e px b>(dx) < oo. 

K |x|>l 

The graphical representation of the functions the Levy measure must 
integrate so that a Levy process has finite moments is given in Figure 7.9. 
The NIG process possesses moments of all order, while the a-stable does not; 
one can already observe in Figure 7.7 that the tails of the Levy measure of 
the a-stable are much heavier than the tails of the NIG. 

Remark 7.4. As can be observed from Propositions 7.1, 7.2 and 7.3, the 
variation of a Levy process depends on the small jumps (and the Brownian 
motion), the moment properties depend on the big jumps, while the activity 
of a Levy process depends on all the jumps of the process. 
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Figure 7.9. A Levy process has first moment if the Levy 
measure integrates \x\ for \x\ > 1 (blue line) and second 
moment if it integrates x 2 for |x| > 1 (orange line). 

8. Some classes of particular interest 

We already know that a Brownian motion, a (compound) Poisson process 
and a Levy jump-diffusion are Levy processes, their Levy-Ito decomposition 
and their characteristic functions. Here, we present some further subclasses 
of Levy processes that are of special interest. 

8.1. Subordinator. A subordinator is an a.s. increasing (in t) Levy process. 
Equivalently, for L to be a subordinator, the triplet must satisfy u{— oo, 0) = 
0, c = 0, j, -^ xu(dx) < oo and 7 = b — f, Q ^ xv{dx) > 0. 
The Levy-Ito decomposition of a subordinator is 



U 



U 



"•/ + / / xn L (ds,dx) 
(0,00) 
and the Levy-Khintchine formula takes the form 

Me iuLt 



•2) 



exp 



lMdx)) 



t(iuj + l (c' 

(0,00) 

Two examples of subordinators are the Poisson and the inverse Gaussian 
process, cf. Figures 8.10 and A. 14. 

8.2. Jumps of finite variation. A Levy process has jumps of finite vari- 
ation if and only if Ji x i <1 \x\v(dx) < 00. In this case, the Levy-Ito decompo- 
sition of L resumes the form 

t 

(8.3) Lt = 7* + VcW t +11 xfi L {ds, dx) 



and the Levy-Khintchine formula takes the form 



U) 



E[e 



iuLt 



exp 



t(iwy 



u 2 c 



+ / (e iux - l)i/(dx)) 
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where 7 is defined similarly to subsection 8.1. 

Moreover, if v([— 1, 1]) < 00, which means that u(M) < 00, then the jumps 
of L correspond to a compound Poisson process. 

8.3. Spectrally one-sided. A Levy processes is called spectrally negative 
if i^(0, 00) = 0. The Levy-Ito decomposition of a spectrally negative Levy 
process has the form 

t t 

(8.5) L t = bt + y/cW t + / / x^ L (ds,dx) + x(fi L - u L )(ds,dx) 

x<-l 0-l<x<0 

and the Levy-Khintchine formula takes the form 

(8.6) JE[e iuLt ] =exp\t(iub-^+ / (e iux - 1 - iul {x> _ 1} )u(dx) 

(-00,0) 

Similarly, a Levy processes is called spectrally positive if — L is spectrally 
negative. 

8.4. Finite first moment. As we have seen already, a Levy process has 
finite first moment if and only if J] , >1 \x\v(dx) < 00. Therefore, we can 
also compensate the big jumps to form a martingale, hence the Levy-Ito 
decomposition of L resumes the form 

t 

(8.7) L t = b't + ^~cW t +11 x( / u L - v L ){ds, dx) 

K 
and the Levy-Khintchine formula takes the form 

(8.8) E[e iuLt ] = exp \t{iub' - ^ + / (e inx - 1 - iux)v{dx)) 

R 

where b' = b + j\ x \ >1 xv{dx). 

Remark 8.1 (Assumption (M)). For the remaining parts we will work 
only with Levy process that have finite first moment. We will refer to them 
as Levy processes that satisfy Assumption (M). For the sake of simplicity, 
we suppress the notation b' and write b instead. 

9. Elements from semimartingale theory 

A semimartingale is a stochastic process X = (Xt)o<t<T which admits 
the decomposition 

(9.1) X = X + M + A 

where Xq is finite and .T-o-measurable, M is a local martingale with Mq = 
and A is a finite variation process with Ao = 0. X is a special semimartingale 
if A is predictable. 

Every special semimartingale X admits the following, so-called, canonical 
decomposition 

(9.2) X = X + B + X c + x*(fi x -u X ). 
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Figure 8.10. Simulated path of a normal inverse Gaussian 
(left) and an inverse Gaussian process. 

Here X c is the continuous martingale part of X and x * (fj, — v ) is the 
•purely discontinuous martingale part of X. fx x is called the random measure 
of jumps of X; it counts the number of jumps of specific size that occur in 
a time interval of specific length. u x is called the compensator of [i x ; for a 
detailed account, we refer to Jacod and Shiryaev (2003, Chapter II). 

Remark 9.1. Note that W * //, for W = W(uj; s, x) and the integer-valued 
measure /i = /j,(ui;dt, dx), t £ [0,T], x £ E, denotes the integral process 



W(uj; t, x)fi(uj; dt, dx). 



E 



Consider a predictable function W : 0, x [0,T] x E 
W * (/j, — v) denotes the stochastic integral 



in Gi oc (n); then 



W (to; t, x) (fj, — v)(u>; dt, dx) . 



E 



Now, recalling the Levy-Ito decomposition (8.7) and comparing it to (9.2), 
we can easily deduce that a Levy process with triplet (6, c, v) which satisfies 
Assumption (M), has the following canonical decomposition 



(9.3) 
where 



and 



E 



L t = bt + y/cWt + / / x(n L - i/ L )(ds,dx), 

b 



/ / xfi L (ds,dx) = ^2 AL 



0<s<t 



x/j, (ds,dx) 



xv (d,s, dx) = t 



Jx.(dx). 
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Therefore, a Levy process that satisfies Assumption (M) is a special semi- 
martingale where the continuous martingale part is a Brownian motion with 
coefficient yfc and the random measure of the jumps is a Poisson random 
measure. The compensator v L of the Poisson random measure \x L is a prod- 
uct measure of the Levy measure with the Lebesgue measure, i.e. v = v%\\ 
one then also writes v (ds, dx) = v{dx)ds. 

We denote the continuous martingale part of L by L c and the purely 
discontinuous martingale part of L by L , i.e. 

(9.4) L r t = ^cW t and Lf = I I x([i L - v L ){ds, dx). 



Remark 9.2. Every Levy process is also a semimartingale; this follows 
easily from (9.1) and the Levy-Ito decomposition of a Levy process. Every 
Levy process with finite first moment (i.e. that satisfies Assumption (M)) 
is also a special semimartingale; conversely, every Levy process that is a 
special semimartingale, has a finite first moment. This is the subject of the 
next result. 

Lemma 9.3. Let L be a Levy process with triplet (b,c,u). The following 
conditions are equivalent 

(1) L is a special semimartingale, 

(2) / R (|x| A |x|>(dx) < oo, 

( 3 ) L\ X \ 1 {\X\>1}V{&X) < OO. 

Proof. From Lemma 2.8 in Kallsen and Shiryaev (2002) we have that, a 
Levy process (semimartingale) is special if and only if the compensator of 
its jump measure satisfies 

(|x| A \x\ 2 )is L (ds,dx) £ V. 
K 

For a fixed t £ R, we get 

t t 

( W AMV(dMx) = //(MAM>(d*)d s 

OR OR 

|2\ 



= t- (\x\ A \x\)u{dx) 

and the last expression is an element of V if and only if 

\x\ A \x\ )v(dx) < oo; 



this settles (1) 44> (2). The equivalence (2) 44> (3) follows from the properties 
of the Levy measure, namely that J\ x \ <x \x\ 2 v(dx) < oo, cf. (7.1). □ 
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10. Martingales and Levy processes 

We give a condition for a Levy process to be a martingale and discuss 
when the exponential of a Levy process is a martingale. 

Proposition 10.1. Let L = (I>t)o<t<T be a Levy process with Levy triplet 
(b,c,v) and assume that E[|Lt|] < oo, i.e. Assumption (M) holds. L is a 
martingale if and only if b = 0. Similarly, L is a submartingale if b > and 
a supermar ting ale if b < 0. 

Proof. The assertion follows immediately from the decomposition of a Levy 
process with finite first moment into a finite variation process, a continuous 
martingale and a pure-jump martingale, cf. equation (9.3). □ 

Proposition 10.2. Let L = (Lt)o<t<T be a Levy process with Levy triplet 
(6, c, v), assume that J, , >1 e ux v{dx) < oo, for net and denote by k the cu- 

mulant of L\, i.e. n{u) = logE[e nLl ]. The process M = (Mt)o<t<T, defined 
via 

R uL t 

is a martingale. 

Proof. Applying Proposition 7.3, we get that E[e" L '] = e tK ( u > < oo, for all 
< t < T. Now, for < s < t, we can re-write M as 

e uL s e u(L t —L s ) e u(Lt—L s ) 

M t = rr u r-r^i = M s n v~rr - 

pSK(U) p(t—S)K(U) p{t — S)K,(U) 

Using the fact that a Levy process has stationary and independent incre- 
ments, we can conclude 

M t 



E 



T s 



MM 



?s 



M e ( t - s ) K ( u ) e ~( t ~ s ) K ( u ) 



L e (t— s)k(w) 

M,. □ 



l s- 



The stochastic exponential 6(L) of a Levy process L = (I*t)o<t<T is the 
solution Z of the stochastic differential equation 

(10.1) dZ t = Z t ^dL t , Z = 1, 
also written as 

(10.2) Z = l + Z_-L, 

where FY means the stochastic integral L F s dY s . The stochastic exponen- 
tial is defined as 



(10.3) £ (L) t = exp U t - \(L c )t) ]J (l + AL 



e -AL 3 



0<s<t 



Remark 10.3. The stochastic exponential of a Levy process that is a mar- 
tingale is a local martingale (cf. Jacod and Shiryaev 2003, Theorem 1.4.61) 
and indeed a (true) martingale when working in a finite time horizon (cf. 
Kallsen 2000, Lemma 4.4). 
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The converse of the stochastic exponential is the stochastic logarithm, 
denoted CogX; for a process X = (X t )o<t<T, the stochastic logarithm is 
the solution of the stochastic differential equation: 

t 
dX 
(10.4) CogX t " 



X S J 
o 



also written as 

(10.5) CogX=^--X. 

Now, if X is a positive process with Xq = 1 we have for Cog X 

(10.6) C 0g X _ log* + -^ ■ m - £ (log (l + ^) - g) 



0<s<- 



for more details see Kallsen and Shiryaev (2002) or Jacod and Shiryaev 
(2003). 

11. Ito's formula 

We state a version of Ito's formula directly for semimartingales, since this 
is the natural framework to work into. 

Lemma 11.1. Let X = (Xt)o<t<T be a real-valued semimartingale and f a 
class C 2 function on M. Then, f(X) is a semimartingale and we have 

t t 

(11.1) f(X t ) = f(X ) + J f'(X s .)dX s + l -f f"(X s _)d(X c ) s 

o o 

+ J2 (f(X s )-f(X s .)-f'(X s .)AX s ), 

0<s<t 

for all t £ [0,T]; alternatively, making use of the random measure of jumps, 
we have 

t t 

(11.2) f(X t ) = f(X ) + j f'(X s .)dX s + l -j f"(X s ^)d(X c ) s 

o o 

t 

+ J j (/(X s _ +x) - f(X 8 -) - f'(X s ^xy x (ds,dx). 

K 

Proof. See Theorem 1.4.57 in Jacod and Shiryaev (2003). □ 

Remark 11.2. An interesting account (and proof) of Ito's formula for Levy 
processes of finite variation can be found in Kyprianou (2006, Chapter 4). 

Lemma 11.3 (Integration by parts). Let X,Y be semimartingales. Then 
XY is also a semimartingale and 

(11.3) XY = I X_dY + I Y-&X + [X,Y\, 
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where the quadratic covariation of X and Y is given by 
(11.4) [X, Y] = (X c , Y c ) + Y, ^X S AY S . 

Proof. See Corollary II. 6. 2 in Protter (2004) and Theorem 1.4.52 in Jacod 
and Shiryaev (2003). □ 

As a simple application of Ito's formula for Levy processes, we will work 
out the dynamics of the stochastic logarithm of a Levy process. 

Let L = (Lt)o<t<T be a Levy process with triplet (6, c, v) and Lq = 1. 
Consider the C 2 function / : M. — ► M. with f(x) = log|x|; then, f'(x) = - 
and f"(x) = — j. Applying Ito's formula to f(L) = log \L\, we get 

t t 

log \L t \ = log \L \ + / — dL s -- j^d(L c ) s 
o o s 

+ Y ( 1 °S\L S \ -log|L s _| - — AL t 



0<s<t s 



t 



^£ogL t = log|L t | + i/^- Y, (lo. 



2 7 Li 

o 0<s<t 



L 



^_ 



AL 



L s - 



Now, making again use of the random measure of jumps of the process L 
and using also that d{L c ) s = &{yfcW) s = cds, we can conclude that 



/ 



c f ds 



CogL t =log\L t \ + - JjT- ( lo S 1 +2 —)^(ds,dx). 



12. GlRSANOV'S THEOREM 

We will describe a special case of Girsanov's theorem for semimartin- 
gales, where a Levy process remains a process with independent increments 
(PII) under the new measure. Here we will restrict ourselves to a finite time 
horizon, i.e. T G [0, oo). 

Let P and P be probability measures defined on the filtered probability 
space (Q,J-, F). Two measures P and P are equivalent, if P(A) = 44> 
P(A) = 0, for all A & J-, and then one writes P ~ P. 

Given two equivalent measures P and P, there exists a unique, positive, 
P-martingale Z = (Zi)o<«r such that Z^ = Ef^lPi], V < t < T. Z is 
called the density process of P with respect to P. 

Conversely, given a measure P and a positive P-martingale Z = (Zt)o<t<T, 
one can define a measure P on (0,P, F) equivalent to P, using the Radon- 
Nikodym derivative E[gp|py] = Zt- 

Theorem 12.1. Let L = (Lt)o<t<T be a Levy process with triplet (b,c,v) 
under P, that satisfies Assumption (M), cf. Remark 8.1. Then, L has the 
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canonical decomposition 

t 

(12.1) L t = bt + ^fcW t + f I ' x(n L -v L )(ds,dx). 

K 

(Al): Assume that P ~ P with density process Z. Then, there exist a 
deterministic process (3 and a measurable non-negative deterministic 
process Y , satisfying 
t 

(12.2) // W y-( S ,x)-l) M dx)d s <oo, 

M 
and 

t 

(c • 0g)ds < oo, 
o 
P-a.s. for < t <T; they are defined by the following formulae: 



(12.3) 



(12.4) 



(Z C ,L C ) = J(c-f3 s -Z s _)ds 



and 



*-"!>(£ 



v\. 



(A2): Conversely, if Z is a positive martingale of the form 

1 



(12.5) 



Z = exp 



p s Vc~dW s - - J ftcds 





+ / J(Y(s..r)- \.)(j, L - v L )(ds.d.r) 




(Y(s, x)-l- ln(Y(s, x)))n L (ds, dx) 



then it defines a probability measure P on (0,, J 7 , F), such that P ~ P. 
(A3): In both cases, we have that W = W — J yfc(3 s ds is a P-Brownian 
motion, v L (ds,dx) = Y(s,x)u L (ds,dx) is the P -compensator of /j, l 
and L has the following canonical decomposition under P: 

t 

(12.6) L t = bt + ^fcW t + I x(fi L -i? L )(ds,dx), 

K 

where 

t t 

(12.7) bt = bt+ / cf3 s ds+ I x(Y(s,x) - l)u L (ds,dx). 

OR 
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Proof. Theorems III. 3. 24, III. 5. 19 and III. 5. 35 in Jacod and Shiryaev (2003) 
yield the result. □ 

Remark 12.2. In (12.4) V = V <g) B(R) is the a-field of predictable sets in 
S7 = S7 x [0,T] xl and M L = \i (uj;dt,dx)P(duj) is the positive measure 
on (fi x [0,T] xl,f®6([0,T])®B(l)) defined by 

(12.8) M* L (W) = E(W*fi L ) T , 

for measurable nonnegative functions W = W(co; t, x) given on 0, x [0, T] x R. 

Z 



Now, the conditional expectation M P L I -y—\V ) is, by definition, the M^-a.s. 
unique "P-measurable function Y with the property 

(12-9) M£(|-t/)=M£(y[/), 

for all nonnegative "P-measurable functions U = U(uj; t, x). 

Remark 12.3. Notice that from condition (12.2) and assumption (M), fol- 
lows that L has finite first moment under P as well, i.e. 

(12.10) E|L t | < oo, for all < t < T. 
Verification follows from Proposition 7.3 and direct calculations. 

Remark 12.4. In general, L is not necessarily a Levy process under the 
measure P; this depends on the tuple (/3, Y). The following cases exist. 

(Gl): if (/3, Y) are deterministic and independent of time, then L re- 
mains a Levy process under P; its triplet is (6, c, Y ■ v). 

(G2): if (P, Y) are deterministic but depend on time, then L becomes 
a process with independent (but not stationary) increments under 
P, often called an additive process. 

(G3): if (/3, Y) are neither deterministic nor independent of time, then 
we just know that L is a semimartingale under P. 

Remark 12.5. Notice that c, the diffusion coefficient, and /j, , the random 
measure of jumps of L, did not change under the change of measure from P 
to P. That happens because c and fi L are path properties of the process and 
do not change under an equivalent change of measure. Intuitively speaking, 
the paths do not change, the probability of certain paths occurring changes. 

Example 12.6. Assume that L is a Levy process with canonical decompo- 
sition (12.1) under P. Assume that P ~ P and the density process is 

t 

(12.11) Z t = exp (3^fcW t +11 ax(fi L - u L )(ds, dx) 

K 



where (3 £ M^o an d a G M. are constants. 

Then, comparing (12.11) with (12.5), we have that the tuple of functions 
that characterize the change of measure is (/?, Y) = (/3, /), where f{x) = e ax . 
Because (/3, f) are deterministic and independent of time, L remains a Levy 
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process under P, its Levy triplet is (6, c, jv) and its canonical decomposition 
is given by equations (12.6) and (12.7). 

Actually, the change of measure of the previous example corresponds to 
the so-called Esscher transformation or exponential tilting. In chapter 3 of 
Kyprianou (2006), one can find a significantly easier proof of Girsanov's 
theorem for Levy processes for the special case of the Esscher transform. 
Here, we reformulate the result of example 12.6 and give a complete proof 
(inspired by Eberlein and Papapantoleon 2005). 

Proposition 12.7. Let L = (L<)o<t<T be a Levy process with canonical 
decomposition (12.1) under P and assume that E[e MZ/t ] < oo for all u £ 
[— p,p], p > 0. Assume that P ~ P with density process Z = (Z t )o<t<T, 
or conversely, assume that P is defined via the Radon-Nikodym derivative 
gp = Zt; here, we have that 



(12.12) 



Z t 



e PL$ e aL? 

E[e^* c ]E[e aL t] 



for (3 £ M. and \a\ < p. Then, L remains a Levy process under P, its Levy 



triplet is (b, c, u), where v = f ■ v for f(x) - 
position is given by the following equations 



and its canonical decom- 



(12.13) 



L t =bt + yfcWt + 



x(fi L - u L )(ds,dx), 



and 
(12.14) 



b + (3c+ / x(e 



l)i/(dx). 



Proof. Firstly, using Proposition 10.2, we can immediately deduce that Z 
is a positive P-martingale; moreover, Zq = 1. Hence, Z serves as a density 
process. 

Secondly, we will show that L has independent and stationary increments 
under P. Using that L has independent and stationary increments under P 
and that Z is a P-martingale, we arrive at the following helpful conclusions: 
for any B £ B(R), F s £ T s and < s < t < T 

(1) l {Lt-L s &B}^- is independent of 1{f 3 } z s and of Z s ; 

(2) E[Z S ] = 1. 

Then, we have that 



P({L t -L s £B}nF s 



E 
E 
E 
E 



1 {L t -L s eB}^{F a }Z t 

Zt 



E [1 W £»] 



1 {L t -L a &B}-^- 

Z i 
l {L t -L 8 eB}TT E[Z s ]E[l {Fa} Z s 

l {L t -L s eB}Zt E[l {Fs }Z s ] 



P({L t -L s £ B})P(F S 
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which yields the independence of the increments. Similarly, regarding the 
stationarity of the increments of L under P, we have that 

P({L t -L s eB})=TE[l {Lt _ LaeB} Z t ] 

= E[l {Lt _ LseB} |-]E[Z s ] 

e a(L- t -L-) e /3(Lf-L d s ) 



E 



E 



l {L t -L s eB} 



h 



^ e a(L<i-Lc) e f3(L?~L^ 
e aL-_ s +f3Lf_ s 



{L! ' LseB1 E[e^- + ^-»] 
= E[l {Lt _ seB} Z t _ s ] 
= P({L t _ s € B}) 

which yields the stationarity of the increments. 

Thirdly, we determine the characteristic function of L under P, which also 
yields the triplet and canonical decomposition. Applying Theorem 25.17 in 
Sato (1999), the moment generating function M^ t of Lt exists for u £ C 
with 3fJit G [— p,p\. We get 



E[e 2Lt ] =TE[e zLt Z t 



E 



e zLt e f3L$ e aLf 

E[e^]E[e aL t 



M e zbt e (z+f3)Lf e (z+a)L?j 

E[eP L t]E[e aL t] 



exp t 



zb + 



(z + (3) 2 c 



+ / ( e ( 2+Q ) x - 1 - (z + a)x)u(dx) 



2 



1 — ax)v(dx) 



exp t 



b + j3c+ / x(e ax -l)i/(dx)) + 



z 2 c 



+ / (e zx - 1 - zx)e ax u{dx) 



exp £ 



- z c 
z6+ — + / [v 



1 — zxW(dx) 



Finally, the statement follows by proving that is(dx) = e ax v{dx) is a Levy 
measure, i.e. / R (l A x 2 )e aa V(dx) < oo. It suffices to note that 



(12.15) 



x 2 e ax is(dx) < C / x 2 z/(dx) < oo, 
|x|<l |x|<l 



where C is a positive constant, because v is a Levy measure; the other part 
follows from the assumptions, since \a\ < p. □ 
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Remark 12.8. Girsanov's theorem is a very powerful tool, widely used in 
mathematical finance. In the second part, it will provide the link between 
the 'real-world' and the 'risk-neutral' measure in a Levy-driven asset price 
model. Other applications of Girsanov's theorem allow to simplify certain 
valuation problems, cf. e.g. Papapantoleon (2007) and references therein. 

13. Construction of Levy processes 

Three popular methods to construct a Levy process are described below. 
(CI): Specifying a Levy triplet; more specifically, whether there ex- 
ists a Brownian component or not and what is the Levy measure. 
Examples of Levy process constructed this way include the stan- 
dard Brownian motion, which has Levy triplet (0, 1, 0) and the Levy 
jump-diffusion, which has Levy triplet (b,a 2 ,XF). 
(C2): Specifying an infinitely divisible random variable as the density 
of the increments at time scale 1 (i.e. L\). Examples of Levy process 
constructed this way include the standard Brownian motion, where 
L\ ~ Normal(0, 1) and the normal inverse Gaussian process, where 
In ~NIG(a,/0,<S,ju). 
(C3): Time-changing Brownian motion with an independent increas- 
ing Levy process. Let W denote the standard Brownian motion; we 
can construct a Levy process by 'replacing' the (calendar) time t 
by an independent increasing Levy process r, therefore L t := W r m, 
< t < T. The process r has the useful - in Finance - interpretation 
as 'business time'. Models constructed this way include the normal 
inverse Gaussian process, where Brownian motion is time-changed 
with the inverse Gaussian process and the variance gamma process, 
where Brownian motion is time-changed with the gamma process. 
Naturally, some processes can be constructed using more than one meth- 
ods. Nevertheless, each method has some distinctive advantages which are 
very useful in applications. The advantages of specifying a triplet (CI) are 
that the characteristic function and the pathwise properties are known and 
allows the construction of a rich variety of models; the drawbacks are that 
parameter estimation and simulation (in the infinite activity case) can be 
quite involved. The second method (C2) allows the easy estimation and sim- 
ulation of the process; on the contrary the structure of the paths might be 
unknown. The method of time-changes (C3) allows for easy simulation, yet 
estimation might be quite difficult. 

14. Simulation of Levy processes 

We shall briefly describe simulation methods for Levy processes. Our at- 
tention is focused on finite activity Levy processes (i.e. Levy jump-diffusions) 
and some special cases of infinite activity Levy processes, namely the nor- 
mal inverse Gaussian and the variance gamma processes. Several speed-up 
methods for the Monte Carlo simulation of Levy processes are presented in 
Webber (2005). 

Here, we do not discuss simulation methods for random variables with 
known density; various algorithms can be found in Devroye (1986), also avail- 
able online at http : //eg . scs . carleton . ca/ luc/rnbookindex . html. 
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14.1. Finite activity. Assume we want to simulate the Levy jump-diffusion 

N t 

L t = bt + aW t + ^Jk 

k=i 

where Nt ~ Poisson(At) and J ~ F(dx). W denotes a standard Brownian 
motion, i.e. Wt ~ Normal(0,t). 

We can simulate a discretized trajectory of the Levy jump-diffusion L at 
fixed time points t\, . . . ,t n as follows: 

• generate a standard normal variate and transform it into a normal 
variate, denoted Gi, with variance oAij, where Ati = U — U—i; 

• generate a Poisson random variate iV with parameter AT; 

• generate iV random variates t& uniformly distributed in [0, T] ; these 
variates correspond to the jump times; 

• simulate the law of jump size J, i.e. simulate random variates Jk 
with law F{dx). 

The discretized trajectory is 

i N 

L U =bU + Y^ G i + Yl 1 {Tk<U} J k- 
i=i fc=i 

14.2. Infinite activity. The variance gamma and the normal inverse Gauss- 
ian process can be easily simulated because they are time-changed Brownian 
motions; we follow Cont and Tankov (2003) closely. A general treatment of 
simulation methods for infinite activity Levy processes can be found in Cont 
and Tankov (2003) and Schoutens (2003). 

Assume we want to simulate a normal inverse Gaussian (NIG) process 
with parameters a, (3, 5, fx; cf. also section 16.5. We can simulate a discretized 
trajectory at fixed time points t\, . . . , t n as follows: 

• simulate n independent inverse Gaussian variables I\ with parame- 
ters (<5Atj) 2 and a 2 — f3 2 , where A£j = U — U-i, i = 1, . . . , n; 

• simulate n i.i.d. standard normal variables Gf, 

• set ALi = fiAU + 13 k + ^/hd. 
The discretized trajectory is 

i 

L ti =^AL fe . 
fc=i 
Assume we want to simulate a variance gamma (VG) process with pa- 
rameters <T, 0, k; we can simulate a discretized trajectory at fixed time points 
t±,...,t n as follows: 

• simulate n independent gamma variables Tj with parameter — i 

• set Ti = nTi] 

• simulate n standard normal variables Gf, 

• set ALi = OTi + o-yT^Gj. 

The discretized trajectory is 

i 

L k = Y, ALk - 
fc=i 
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Part 2. Applications in Finance 

15. Asset price model 

We describe an asset price model driven by a Levy process, both under 
the 'real' and under the 'risk-neutral' measure. Then, we present an informal 
account of market incompleteness. 

15.1. Real-world measure. Under the real-world measure, we model the 
asset price process as the exponential of a Levy process, that is 

(15.1) S t = S expL t , 0<t<T, 

where, L is the Levy process whose infinitely divisible distribution has been 
estimated from the data set available for the particular asset. Hence, the log- 
returns of the model have independent and stationary increments, which are 
distributed - along time intervals of specific length, e.g. 1 - according to an 

infinitely divisible distribution C(X), i.e. L\ = X. 

Naturally, the path properties of the process L carry over to S; if, for 
example, L is a pure-jump Levy process, then S is also a pure-jump process. 
This fact allows us to capture, up to a certain extent, the microstructure of 
price fluctuations, even on an intraday time scale. 

An application of Ito's formula yields that S = (St)o<t<T is the solution 
of the stochastic differential equation 

(15.2) dS t = S t - (dL t + |dt + f(e x - 1 - x)fi L (dt, da;)) . 

K 

We could also specify S by replacing the Brownian motion in the Black- 
Scholes SDE by a Levy process, i.e. via 

(15.3) dSt = S t -dL t , 
whose solution is the stochastic exponential 

(15.4) S t = S £(L t ). 

The second approach is unfavorable for financial applications, because (a) 
the asset price can take negative values, unless jumps are restricted to be 
larger than — 1, i.e. supp(^) C [— 1, oo), and (6) the distribution of log-returns 
is not known. Of course, in the special case of the Black-Scholes model the 
two approaches coincide. 

Remark 15.1. The two modeling approaches are nevertheless closely re- 
lated and, in some sense, complementary of each other. One approach is 
suitable for studying the distributional properties of the price process and 
the other for investigating the martingale properties. For the connection be- 
tween the natural and stochastic exponential for Levy processes, we refer to 
Lemma A.8 in Goll and Kallsen (2000). 

The fact that the price process is driven by a Levy process, makes the 
market, in general, incomplete; the only exceptions are the markets driven 
by the Normal (Black-Scholes model) and Poisson distributions. Therefore, 
there exists a large set of equivalent martingale measures, i.e. candidate 
measures for risk-neutral valuation. 
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Eberlein and Jacod (1997) provide a thorough analysis and characteriza- 
tion of the set of equivalent martingale measures for Levy-driven models. 
Moreover, they prove that the range of option prices for a convex payoff 
function, e.g. a call option, under all possible equivalent martingale mea- 
sures spans the whole no-arbitrage interval, e.g. [(So — Ke~ rT ) + , Sq] for a 
European call option with strike K. Selivanov (2005) discusses the existence 
and uniqueness of martingale measures for exponential Levy models in finite 
and infinite time horizon and for various specifications of the no-arbitrage 
condition. 

The Levy market can be completed using particular assets, such as mo- 
ment derivatives (e.g. variance swaps), and then there exists a unique equiv- 
alent martingale measure; see Corcuera, Nualart, and Schoutens (2005a, 
2005b). For example, if an asset is driven by a Levy jump-diffusion 

N t 

(15.5) L t = bt + J~cW t + J2 J k 

fe=i 

where Jt = ctMk, then the market can be completed using only variance 
swaps on this asset; this example will be revisited in section 15.3. 

15.2. Risk-neutral measure. Under the risk neutral measure, denoted by 
P, we model the asset price process as an exponential Levy process 

(15.6) S t = SoexpL t 

where the Levy process L has the triplet (b, c, v) and satisfies Assumptions 
(M) (cf. Remark 8.1) and (EM) (see below). 
The process L has the canonical decomposition 

t 

(15.7) Lt = It + Vm t + I f x(n L - u L )(ds, dx) 

R 

where W is a P-Brownian motion and v is the P-compensator of the jump 
measure \x L . 

Because we have assumed that P is a risk neutral measure, the asset price 
has mean rate of return \i = r — 5 and the discounted and re-invested process 
( e ( r ~ '")* St)o<t<T j is a martingale under P. Here r > is the (domestic) risk- 
free interest rate, 5 > the continuous dividend yield (or foreign interest 
rate) of the asset. Therefore, the drift term b takes the form 

(15.8) b = r-b- C -- f(e x - 1 - x)v(dx); 

R 

see Eberlein, Papapantoleon, and Shiryaev (2008) and Papapantoleon (2007) 
for all the details. 

Assumption (EM). We assume that the Levy process L has finite first 
exponential moment, i.e. 

(15.9) E[e Lt ] < oo. 
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There are various ways to choose the martingale measure such that it 
is equivalent to the real-world measure. We refer to Goll and Ruschendorf 
(2001) for a unified exposition - in terms of /-divergences - of the different 
methods for selecting an equivalent martingale measure (EMM). Note that, 
some of the proposed methods to choose an EMM preserve the Levy property 
of log-returns; examples are the Esscher transformation and the minimal 
entropy martingale measure (cf. Esche and Schweizer 2005). 

The market practice is to consider the choice of the martingale measure 
as the result of a calibration to market data of vanilla options. Hakala and 
Wystup (2002) describe the calibration procedure in detail. Cont and Tankov 
(2004, 2006) and Belomestny and Reifi (2005) present numerically stable 
calibration methods for Levy driven models. 

15.3. On market incompleteness. In order to gain a better understand- 
ing of why the market is incomplete, let us make the following observation. 
Assume that the price process of a financial asset is modeled as an exponen- 
tial Levy process under both the real and the risk-neutral measure. Assume 
that these measures, denoted P and P, are equivalent and denote the triplet 
of the Levy process under P and P by (6, c, v) and (b, c, v) respectively. 

Now, applying Girsanov's theorem we get that these triplets are related 
via c = c, v = Y ■ v and 

(15.10) b = b + cf3 + x(Y -l)*u, 

where ((3, Y) is the tuple of functions related to the density process. On the 
other hand, from the martingale condition we get that 

(15.11) b = r---(e x -l-x)*i?. 

Equating (15.10) and (15.11) and using c = c and v = Y ■ D, we have that 
= b + c/3 + x(Y -l)*v-r+- + (e x -l-x)*v 

(15.12) «• = b - r + c(/3 + -) + ((e x - 1)Y - x) * v; 

therefore, we have one equation but two unknown parameters, (3 and Y 
stemming from the change of measure. Every solution tuple (/?, Y) of equa- 
tion (15.12) corresponds to a different equivalent martingale measure, which 
explains why the market is not complete. The tuple (/3, 1") could also be 
termed the tuple of 'market price of risk '. 

Example 15.2 (Black-Scholes model). Let us consider the Black-Scholes 
model, where the driving process is a Brownian motion with drift, i.e. Lt = 
bt + \/cWt. Then, equation (15.12) has a unique solution, namely 

r — b 1 

(15.13) £=___, 

the martingale measure is unique and the market is complete. We can also 
easily check that plugging (15.13) into (15.10), we recover the martingale 
condition (15.11). 
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Remark 15.3. The quantity (3 in (15.13) is nothing else than the so-called 
market price of risk. The difference from the quantity often encountered in 
textbooks, i.e. Tj= ^, stems from the fact that we model using the natural 
instead of the stochastic exponential, i.e. using SDE (15.2) and not (15.3). 

Example 15.4 (Poisson model). Let us consider the Poisson model, where 
the driving motion is a Poisson process with intensity A > and jump size 
a, i.e. Lt = bt + aNt and v(dx) = Xli a \(dx). Then, equation (15.12) has a 
unique solution for Y, which is 

= b - r + ((e x - 1)Y - x) * Al {o} (dx) 

& = b - r + ((e a - l)Y - a)X 

(15.14) ^Y = - — — ; 

y ' (e a - 1)A 

therefore the martingale measure is unique and the market is complete. By 
the analogy to the Black-Scholes case, we could call the quantity Y in (15.14) 
the market price of jump risk. 

Moreover, we can also check that plugging (15.14) into (15.10), we recover 
the martingale condition (15.11); indeed, we have that 

b = b + a\(Y-l) 

= b + aX{Y - 1) + {e a - 1)YX - {e a - 1)YX 

= r-{e a -l- a)X, 

where we have used (15.14) and that P = Y ■ is, which in the current frame- 
work translates to A = YX. 

Example 15.5 (A simple incomplete model). Assume that the driving 
process consists of a drift, a Brownian motion and a Poisson process, i.e. 
Lt = bt + \fcWt + aNt, as in examples 15.2 and 15.4. Based on (15.13) and 

(15.14) we postulate that the solutions of equation (15.12) are of the form 

(15.15) P e = e r -^- 1 and Y £ = <± - f r - % + " A 
V ' ' c 2 (e Q -l)A 

for any e £ (0, 1). One can easily verify that f3 £ and Y £ satisfy (15.12). But 
then, to any e G (0, 1) corresponds an equivalent martingale measure and 
we can easily conclude that this simple market is incomplete. 

16. Popular models 

In this section, we review some popular models in the mathematical fi- 
nance literature from the point of view of Levy processes. We describe their 
Levy triplets and characteristic functions and provide, whenever possible, 
their - infinitely divisible - laws. 

16.1. Black Scholes. The most famous asset price model based on a Levy 
process is that of Samuelson (1965), Black and Scholes (1973) and Merton 
(1973). The log-returns are normally distributed with mean \x and variance 
a 2 , i.e. L\ ~ Normal(/i, a 2 ) and the density is 

, / n 1 r (x- fi) 2 ] 

OV27T L 2fT J 



a 2 u 2 



34 ANTONIS PAPAPANTOLEON 

The characteristic function is 

fLi ( u ) = exp ifiu 
the first and second moments are 

E[Li] = /z, Var[Lx] = a 2 , 
while the skewness and kurtosis are 

skew[Li] = 0, kurt[Li] = 3. 
The canonical decomposition of L is 

Lt = nt + <rW t 
and the Levy triplet is (//, a 2 , 0). 

16.2. Merton. Merton (1976) was one of the first to use a discontinuous 
price process to model asset returns. The canonical decomposition of the 
driving process is 

Nt 

Lt = fd + aW t + Y, J k 

fc=i 

where Jk ~ Normal(// j,0j), k = 1, ..., hence the distribution of the jump 
size has density 



fj(x) = -±= 

The characteristic function of L\ is 



exp 



(x-fijf 



2a 2 



(p Ll (u) = exp 



a 2 u 2 



IflU 



+ \U 1 ^j u ~ u2 j' u - 2 / 2 — i) 



and the Levy triplet is (//, a 2 , A x fj). 

The density of L\ is not known in closed form, while the first two moments 
are 

Var[Li] = a 2 + A/ij + Xa 2 



E[Li] = n + A//j 



and 



16.3. Kou. Kou (2002) proposed a jump-diffusion model similar to Mer- 
ton's, where the jump size is double-exponentially distributed. Therefore, 
the canonical decomposition of the driving process is 

N t 

L t = fit + aW t + Y, J k 

fc=i 

where Jk ~ DbExpo(p, d\, O2), k = 1, ..., hence the distribution of the jump 
size has density 

fj[x) =p9 1 e~ e ^l {x<0} + (1 -p)d 2 e 02X l {x>o} . 

The characteristic function of L\ is 

P0l (l-p)02 



V?Li («) = exp 



aV 



IjXU 



+ A 



71 — IU 



02 + 2W 



1 



and the Levy triplet is (//, <r 2 , A x /j). 
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The density of L\ is not known in closed form, while the first two moments 



arc 



V\T 1 , X P K 1 -?) i ,,„ Ir , 2 . Ap A(i-p) 

E[Li] = ix + and Var[Li] = a + -^ + — ^ ■ 

16.4. Generalized Hyperbolic. The generalized hyperbolic model was 
introduced by Eberlein and Prause (2002) following the seminal work on 
the hyperbolic model by Eberlein and Keller (1995). The class of hyper- 
bolic distributions was invented by O. E. Barndorff-Nielsen in relation to 
the so-called 'sand project' (cf. Barndorff-Nielsen 1977). The increments of 
time length 1 follow a generalized hyperbolic distribution with parameters 
a, (3, 5, fi, A, i.e. L\ ~ GH(a, (3, 5, fx, A) and the density is 

f GH (x) = c(A, a, (3, 5)(5 2 + (x - li) 2 ) (X ^ )/2 

x K x _i (a\f5 2 + (x — /j,) 2 ) exp (/3(x — //)), 

where 



c(A, a, (3, 5) 



a 



2 



/3 2 ) A /2 



2ira x -*K x (6y/a 2 -p 2 ) 



and K\ denotes the Bessel function of the third kind with index A (cf. 
Abramowitz and Stegun 1968). Parameter a > determines the shape, 
< |/3 1 < a determines the skewness, (i£l the location and 5 > is a 
scaling parameter. The last parameter, A € R affects the heaviness of the 
tails and allows us to navigate through different subclasses. For example, for 
A = 1 we get the hyperbolic distribution and for X = —^ we get the normal 
inverse Gaussian (NIG). 

The characteristic function of the GH distribution is 

while the first and second moments are 

(36 2 K X+1 (Q 

and 

VarfLx] = S l ^±±tt} + ^ ^A +2 (C) ^A+i(0 



C K X (Q C V K X (Q kko 

where (, = 5\/a 2 — [3 2 . 

The canonical decomposition of a Levy process driven by a generalized 
hyperbolic distribution (i.e. L\ ~ GH) is 

t 
L t = tE[Li] + j I x(n L - v GH )(ds,dx) 

R 
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and the Levy triplet is (E[Li],0, v ). The Levy measure of the GH distri- 
bution has the following form 

GHia n ePx ( f exp(- v / 2y + a 2 |x|) _ a]x] \ 

here J^ and 1^ denote the Bessel functions of the first and second kind with 
index A. We refer to Raible (2000, section 2.4.1) for a fine analysis of this 
Levy measure. 

The GH distribution contains as special or limiting cases several known 
distributions, including the normal, exponential, gamma, variance gamma, 
hyperbolic and normal inverse Gaussian distributions; we refer to Eberlein 
and v. Hammerstein (2004) for an exhaustive survey. 

16.5. Normal Inverse Gaussian. The normal inverse Gaussian distribu- 
tion is a special case of the GH for A = — ^ it was introduced to finance in 
Barndorff-Nielsen (1997). The density is 



a KAad^l + i^) 2 

Inig(x) = - exp (6y/a 2 - (3 2 + (3(x - n)) — . 

while the characteristic function has the simplified form 

<PMG(u) = e^ g ^gEg L. 

exp((5y / a 2 — (j3 + iu) 2 ) 

The first and second moments of the NIG distribution are 
E[L 1 ] = n+— |L= and Var^] = -=£== + ^ 



yja 2 - (3 2 ^/a 2 - (3 2 (^/a 2 - (3 2 ) 3 

and similarly to the GH, the canonical decomposition is 

t 
L t = tE[Li] + I I x(fi L - u NIG )(ds,dx), 

K 

where now the Levy measure has the simplified form 
v NIG (dx) = e^^KUalxttdx. 

TT\X\ 

The NIG is the only subclass of the GH that is closed under convolution, 
i.e. if X ~ NIG(a, /?, S±, /Ui) and Y ~ NIG(a, (3, 62, ^2) and X is independent 
of Y, then 

X + Y ~ NIG(a, 0, <Ji + <5 2 , fa + ix 2 ). 

Therefore, if we estimate the returns distribution at some time scale, then 
we know it - in closed form - for all time scales. 
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16.6. CGMY. The CGMY Levy process was introduced by Carr, Geman, 
Madan, and Yor (2002); another name for this process is (generalized) tem- 
pered stable process (see e.g. Cont and Tankov 2003). The characteristic 
function of L t , t G [0, T] is 

ip Lt (u) = exp (tCT(-Y) [{M - iu) Y + (G + iu) Y - M Y - G Y 



The Levy measure of this process admits the representation 

-Mx Gx 

a .I+y 1 {x>0}dx + C|^ 



V CGMY 



e e 

( dx ) = C ZaW 1 {»>o}da; + C—^l^oydx, 



where C > 0, G > 0, M > 0, and Y < 2. The CGMY process is a pure jump 
Levy process with canonical decomposition 



L t = tE[Li] + / / x(/i L - u^ MY )(ds,dx), 



and Levy triplet (E[Li],0, u CGMY ), while the density is not known in closed 
form. 

The CGMY processes are closely related to stable processes; in fact, the 
Levy measure of the CGMY process coincides with the Levy measure of the 
stable process with index a G (0, 2) (cf. Samorodnitsky and Taqqu 1994, Def. 
1.1.6), but with the additional exponential factors; hence the name tempered 
stable processes. Due to the exponential tempering of the Levy measure, 
the CGMY distribution has finite moments of all orders. Again, the class of 
CGMY distributions contains several other distributions as subclasses, for 
example the variance gamma distribution (Madan and Seneta 1990) and the 
bilateral gamma distribution (Kuchlcr and Tappe 2008). 

16.7. Meixner. The Meixner process was introduced by Schoutens and 
Teugels (1998), see also Schoutens (2002). Let L = (Lt)o<t<T be a Meixner 
process with Law (H\\P) = Meixner ( a, /3, 5), a > 0, —it < (3 < n, 5 > 0, 
then the density is 



8\ 2S 

2cos 2j (j3x 



/M ™ (X) = 2avT{25) eXP V o 
The characteristic function Lt, t £ [0,T] is 



a 



l \ I C ° S 2 

<PLt{u) ~ 



261 



cosh 



au—i/3 



2 / 

and the Levy measure of the Meixner process admits the representation 

5 exp ( !-x 
i/ Meixncr (d^ - 



xsinh( 



TT.r > 
a ' 
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The Meixner process is a pure jump Levy process with canonical decompo- 
sition 

t 

, L , .Mcixncr \ 



L t = ffi[Li] + / / x(fi L - z/ Viclxncr )(ds,chr), 
K 

and Levy triplet (E[Li],0, z/ Mcixner ). 

17. Pricing European options 

The aim of this section is to review the three predominant methods 
for pricing European options on assets driven by general Levy processes. 
Namely, we review transform methods, partial integro-differential equation 
(PIDE) methods and Monte Carlo methods. Of course, all these methods 
can be used - under certain modifications - when considering more general 
driving processes as well. 

The setting is as follows: we consider an asset S = (St)o<t<T modeled as 
an exponential Levy process, i.e. 

(17.1) S t = S expL t , 0<t<T, 

where L = (Xt)o<t<T has the Levy triplet (6, c, v). We assume that the 
asset is modeled directly under a martingale measure, cf. section 15.2, hence 
the martingale restriction on the drift term b is in force. For simplicity, we 
assume that r > and 6 = throughout this section. 

We aim to derive the price of a European option on the asset S with 
payoff function g maturing at time T, i.e. the payoff of the option is g(Sx)- 

17.1. Transform methods. The simpler, faster and most common method 
for pricing European options on assets driven by Levy processes is to de- 
rive an integral representation for the option price using Fourier or Laplace 
transforms. This blends perfectly with Levy processes, since the represen- 
tation involves the characteristic function of the random variables, which is 
explicitly provided by the Levy-Khint chine formula. The resulting integral 
can be computed numerically very easily and fast. The main drawback of 
this method is that exotic derivatives cannot be handled so easily. 

Several authors have derived valuation formulae using Fourier or Laplace 
transforms, see e.g. Carr and Madan (1999), Borovkov and Novikov (2002) 
and Eberlein, Glau, and Papapantoleon (2008). Here, we review the method 
developed by S. Raible (cf. Raible 2000, Chapter 3). 

Assume that the following conditions regarding the driving process of the 
asset and the payoff function are in force. 

(Tl): Assume that (fL T (z), the extended characteristic function of Lt, 
exists for all z € C with Qz £ I\ D [0, 1]. 

(T2): Assume that Pl t , the distribution of Lt, is absolutely continu- 
ous w.r.t. the Lebesgue measure A with density p. 

(T3): Consider an integrable, European-style, payoff function q(St)- 

(T4): Assume that x ^ e~ \g(e~ x )\ is bounded and integrable for 
all R £ I 2 C R. 

(T5): Assume that h (~l h + 0- 
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Furthermore, let £h(z) denote the bilateral Laplace transform of a func- 
tion h at z € C, i.e. let 



£ h (z) := / e~ zx h(x)dx. 

11 
According to arbitrage pricing, the value of an option is equal to its dis- 
counted expected payoff under the risk-neutral measure P. Hence, we get 

C T (S, K) = e- rT n9(S T )] = e~ rT Jg(S T )dP 

h 

J g(S e x )dP LT (x) = e~ rT j g(S e x )p(x)dx 



e~ rT 



because Pl t is absolutely continuous with respect to the Lebesgue measure. 
Define the function ir{x) = g(e~ x ) and let £ = —log So, then 

(17.2) C T (S, K) = e~ rT i vr(C - x)p{x)dx = e^vr * p)(Q =■ C 

'n 
which is a convolution of it with p at the point (, multiplied by the discount 
factor. 

The idea now is to apply a Laplace transform on both sides of (17.2) and 
take advantage of the fact that the Laplace transform of a convolution equals 
the product of the Laplace transforms of the factors. The resulting Laplace 
transforms are easier to calculate analytically. Finally, we can invert the 
Laplace transforms to recover the option value. 

Applying Laplace transforms on both sides of (17.2) for C B z = R + 
iu, R £ I\ n I2, u £ R, we get that 

£c(z) = e~ rT fe- zx (TT * p){x)dx 

R 

= e~ rT / e~ zx n(x)dx / e~ zx p(x)dx 

R R 

= e- rT £ n (z)2 p (z). 
Now, inverting this Laplace transform yields the option value, i.e. 

R+ioo 

°t(S,K) = ^- J e< z £ c (z)dz 

R—ioo 

— I ' S R+iu ^ c {R + iu)du 
2vr J 

f e iC "e- rT £^(i? + iu)£ p (R + iu)du 
/ e iC "£ 7r ( J R + iu)cp LT (iR - u)du. 



R—ioo 

ft 

~2tt 

-rT+CR 



2tt 
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Here, 2,^ is the Laplace transform of the modified payoff function tt(x) = 
g(e~ x ) and <pl t is provided directly from the Levy-Khintchine formula. 
Below, we describe two important examples of payoff functions and their 
Laplace transforms. 

Example 17.1 (Call and put option). A European call option pays off 
giSx) = (St — K) + , for some strike price K. The Laplace transform of its 
modified payoff function n is 

K l+z 

(17 - 3) £ ' W = i^TT) 

for z G C with Viz = R G I 2 = (-oo, -1). 

Similarly, for a European put option that pays off g(Sx) = (K — St) + , 
the Laplace transform of its modified payoff function it is given by (17.3) for 
z G C with $tz = Re h = (0, oo). 

Example 17.2 (Digital option). A European digital call option pays off 
g(Sx) = 1ss t >k\- The Laplace transform of its modified payoff function ir 
is 

(17.4) £ w (z) = ~ 

for z G C with $tz = R e I 2 = (-oo,0). 

Similarly, for a European digital put option that pays offp(S'r) = lig T< x\, 
the Laplace transform of its modified payoff function tt is 

K z 

(17.5) £*(*) = — 

for z G C with $z = R G I 2 = (0, oo). 

17.2. PIDE methods. An alternative to transform methods for pricing 
options is to derive and then solve numerically the partial integro-differential 
equation (PIDE) that the option price satisfies. Note that in their seminal 
paper Black and Scholes derive such a PDE for the price of a European 
option. The advantage of PIDE methods is that complex and exotic payoffs 
can be treated easily; the limitations are the slower speed in comparison 
to transform methods and the computational complexity when handling 
options on several assets. 

Here, we derive the PIDE corresponding to the price of a European option 
in a Levy-driven asset, using martingale techniques; of course, we could 
derive the same PIDE by constructing a self-financing portfolio. 

Let us denote by G(St, t) the time-t price of a European option with payoff 
function g on the asset S; the price is given by 

(17.6) G(S t , t) = e- r ( T -<)E[<?(DY)] =:V t , 0<t< T. 

By arbitrage theory, we know that the discounted option price process must 
be a martingale under a martingale measure. Therefore, any decomposition 
of the price process as 

(17.7) e~ rt V t = V + M t + A t , 

where M G Mi oc and A G A\ oc , must satisfy A t = for all t G [0, T]. This 
condition yields the desired PIDE. 
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Now, for notational but also computational convenience, we work with 
the driving process L and not the asset price process S, hence we derive a 
PIDE involving f(Lf,t) := G(Sf,t), or in other words 

(17.8) f(L t , t) = e-^ T -^n9(S e L n] = V t , < t < T. 

Let us denote by dif the derivative of / with respect to the i-th argument, 
dff the second derivative of / with respect to the i-th argument, and so on. 
Assume that / G C 2,1 (IR x [0, T]), i.e. it is twice continuously differentiable 
in the first argument and once continuously differentiable in the second 
argument. An application of Ito's formula yields: 

d(e- H V t ) = d(e- rt f(L t _,t)) 

= -re- rt f(L t „ ,t)dt + e- rt d 2 f(L t ^ , t)dt 

+ e- rt d 1 f(L t _,t)dL t + h- rt d 2 1 f(L t „,t)dm) 

+ e~ rt J (f(L t _+z,t) - f(L t _,t) - &/(£*_, i)z)M L (d*,d*) 



e -rt J _ r/(L< _ ) t)dt + d 2 f( Lt _ , t )dt + Si/(Lt_ , t)bdt 

+ d 1 f(L t _,t)^dW t + fd 1 f(L t ^,t)z(fi L -i^ L )(dz,dt) 



+ l -dlf{L t _,t)cdt 



+ J[f{L t _ +z,t)- f{Lt-,t) - d!/(L t _ , t)z) (,i L - u L )(dz, dt) 
+ f(f(L t _ +z,t)- f(Lt- , t) - d!f(L t _ , t)z) u(dz)dt) . 



Now, the stochastic differential of the bounded variation part of the option 
price process is 

e- ri {-r/(L t _,t)+a 2 /(L t _,t) + 9 1 /(L t _,t)6+^ 1 2 /(L t _,t)c 



+ J[f(L t _+z,t)-f(L t -,t)-d 1 f(Lt_,t)z)v(dz)>, 

K 

while the remaining parts constitute of the local martingale part. 

As was already mentioned, the bounded variation part vanishes identi- 
cally. Hence, the price of the option satisfies the partial integro-differential 
equation 

(17.9) = -rf(x, t) + 8 2 f(x, t) + 6>!/(x, t)b + |a 2 /(x, t) 

+ f[f{x + z,t)- f(x,t) - d 1 f(x,t) z y(dz), 
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for all (x,t) Eix (0, T), subject to the terminal condition 
(17.10) f(x,T)=g(e x 



■+x\ 



Remark 17.3. Using the martingale condition (15.8) to make the drift term 
explicit, we derive an equivalent formulation of the PIDE: 



C ^ -' - C ^2 



= -rf(x, t) + d 2 f(x, t) + (r- -)0i/(z, t) + -d(f(x, t) 



+ 



(j(x + z,t) - f(x,t) - (e z - l)d 1 f(x,t)y(dz), 



for all (x,t) eix (0,T), subject to the terminal condition 

f(x,T) = g(e x ). 

Remark 17.4. Numerical methods for solving the above partial integro- 
differential equations can be found, for example, in Matache et al. (2004), 
in Matache et al. (2005) and in Cont and Tankov (2003, Chapter 12). 

17.3. Monte Carlo methods. Another method for pricing options is to 
use a Monte Carlo simulation. The main advantage of this method is that 
complex and exotic derivatives can be treated easily - which is very impor- 
tant in applications, since little is known about functionals of Levy processes. 
Moreover, options on several assets can also be handled easily using Monte 
Carlo simulations. The main drawback of Monte Carlo methods is the slow 
computational speed. 

We briefly sketch the pricing of a European call option on a Levy driven 
asset. The payoff of the call option with strike K at the time of maturity T 
is g(ST) = (St — K) + and the price is provided by the discounted expected 
payoff under a risk-neutral measure, i.e. 

C T (S,K) = e- rT TE[(S T - K) + }. 

The crux of pricing European options with Monte Carlo methods is to sim- 
ulate the terminal value of asset price St = 5*0 exp Lt - see section 14 for 
simulation methods for Levy processes. Let Sr fe for k = 1,...,N denote 
the simulated values; then, the option price Ct(S, K) is estimated by the 
average of the prices for the simulated asset values, that is 

N 

C T (S,K)=e- rT J2(ST k -K)+, 
fc=i 
and by the Law of Large Numbers we have that 

C T (S, K) -► C T (S, K) as N -► oo. 
18. Empirical evidence 

Levy processes provide a framework that can easily capture the empiri- 
cal observations both under the "real world" and under the "risk-neutral" 
measure. We provide here some indicative examples. 

Under the "real world" measure, Levy processes are generated by dis- 
tributions that are flexible enough to capture the observed fat-tailed and 
skewed (leptokurtic) behavior of asset returns. One such class of distribu- 
tions is the class of generalized hyperbolic distributions (cf. section 16.4). In 
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Figure 18.11. Densities of hyperbolic (red), NIG (blue) and 
hyperboloid distributions (left). Comparison of the GH (red) 
and Normal distributions (with equal mean and variance). 




Figure 18.12. Empirical distribution and Q-Q plot of 
EUR/USD daily log-returns with fitted GH (red). 



Figure 18.11, various densities of generalized hyperbolic distributions and a 
comparison of the generalized hyperbolic and normal density are plotted. 

A typical example of the behavior of asset returns can be seen in Figures 
1.2 and 18.12. The fitted normal distribution has lower peak, fatter flanks 
and lighter tails than the empirical distribution; this means that, in real- 
ity, tiny and large price movements occur more frequently, and small and 
medium size movements occur less frequently, than predicted by the nor- 
mal distribution. On the other hand, the generalized hyperbolic distribution 
gives a very good statistical fit of the empirical distribution; this is further 
verified by the corresponding Q-Q plot. 

Under the "risk-neutral" measure, the flexibility of the generating dis- 
tributions allows the implied volatility smiles produced by a Levy model 
to accurately capture the shape of the implied volatility smiles observed in 
the market. A typical volatility surface can be seen in Figure 1.3. Figure 
18.13 exhibits the volatility smile of market data (EUR/USD) and the cal- 
ibrated implied volatility smile produced by the NIG distribution; clearly, 
the resulting smile fits the data particularly well. 



11 
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Figure 18.13. Implied volatilities of EUR/USD options and 
calibrated NIG smile. 



Appendix A. Poisson random variables and processes 

Definition A.l. Let X be a Poisson distributed random variable with pa- 
rameter A € M^o- Then, for n £ N the probability distribution is 

A n 



P(X 



11) 



-X- 



IV. 



and the first two centered moments are 

E[X} = X and 



VarLY] = A. 



Definition A. 2. A cadlag, adapted stochastic process N = (Nt)o<t<T with 
N t : SI x M^o ->NU {0} is called a Poisson process if 

(1) N = 0, 

(2) N t - N s is independent of T s for any < s < t < T, 

(3) Nt — N s is Poisson distributed with parameter \(t — s) for any < 
s<t<T. 

Then, A > is called the intensity of the Poisson process. 

Definition A. 3. Let N be a Poisson process with parameter A. We shall 
call the process N = (iVi)o<t<T with Nt '■ Q X M^o — * R where 

(A.l) N t := N t - At 

a compensated Poisson process. 

A simulated path of a Poisson and a compensated Poisson process can be 
seen in Figure A. 14. 

Proposition A. 4. The compensated Poisson process defined by (A.l) is a 
martingale. 

Proof. We have that 

(1) the process N is adapted to the filtration because N is adapted (by 
definition); 

(2) E[|iVt|] < oo because TE,[\N t \] < oo, for all < t < T; 
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Figure A. 14. Plots of the Poisson (left) and compensated 
Poisson process. 



(3) finally, let < s < t < T, then 

E\N t \F a ] = E[N t - \t\F a ] 

= E[iV s - (N t - N S )\F S ] - X(s -(t- s)) 
= N S - E[(JV t - N S )\F S ] - \(s -(t- s)) 
= N S -Xs 



□ 



Remark A. 5. The characteristic functions of the Poisson and compensated 
Poisson random variables are respectively 



and 



E[e iuJVt ] = exp [Xt(e iu - 1)] 
TE[e iuNt ] = exp [\t(e iu - 1 - iu)] . 



Appendix B. Compound Poisson random variables 

Let N be a Poisson distributed random variable with parameter A > 
and J = (Jfc)fc>i an i-i-d. sequence of random variables with law F. Then, 
by conditioning on the number of jumps and using independence, we have 
that the characteristic function of a compound Poisson distributed random 
variable is 



E 



n 



P(N = n) 



n>0 



5~)E[e fa £J= 



=1 Jk 



n>0 



A' 



111 




e lux F(dx) | e" v 



in 



exp A / (e* 



l)F(dx) 
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Appendix C. Notation 

= equality in law, C{X) law of the random variable X 
a Ab = min{a, b}, a V b = max{a, 6} 
C 3 z = a + i(3, with a, (3 £ M; then 3ftz = a and Qz = (3 
1a denotes the indicator of the generic event A, i.e. 

i A / T \ - J !. if are A, 

Classes: 

A'tioc local martingales 

^ll oc processes of locally bounded variation 

V processes of finite variation 

Gioc(m) functions integrable wrt the compensated random measure [i — v 



Appendix D. Datasets 

The EUR/USD implied volatility data are from 5 November 2001. The 
spot price was 0.93, the domestic rate (USD) 5% and the foreign rate (EUR) 
4%. The data are available at 

http://www.mathfinance.de/FF/sampleinputdata.txt. 

The USD/JPY, EUR/USD and GBP/USD foreign exchange time se- 
ries correspond to noon buying rates (dates: 22/10/1997 — 22/10/2004, 
4/1/99 - 3/3/2005 and 1/5/2002 - 3/3/2005 respectively). The data can 
be downloaded from 

http : //www . newyorkf ed . org/markets/f oreignex . html. 



Appendix E. Paul Levy 

Processes with independent and stationary increments are named Levy 
processes after the French mathematician Paul Levy (1886-1971), who made 
the connection with infinitely divisible laws, characterized their distributions 
(Levy-Khintchine formula) and described their path structure (Levy-Ito de- 
composition). Paul Levy is one of the founding fathers of the theory of 
stochastic processes and made major contributions to the field of probabil- 
ity theory. Among others, Paul Levy contributed to the study of Gaussian 
variables and processes, the law of large numbers, the central limit theorem, 
stable laws, infinitely divisible laws and pioneered the study of processes 
with independent and stationary increments. 

More information about Paul Levy and his scientific work, can be found 
at the websites 

http : //www . cmap . polytechnique . f r/ rama/levy . html 
and 

http : //www . annales . org/archives/x/paullevy . html 
(in French). 



introduction to levy processes 47 

Acknowledgments 

A large part of these notes was written while I was a Ph.D. student at 
the University of Freiburg; I am grateful to Ernst Eberlein for various inter- 
esting and illuminating discussions, and for the opportunity to present this 
material at several occasions. I am grateful to the various readers for their 
comments, corrections and suggestions. Financial support from the Deutsche 
Forschungsgemeinschaft (DFG, Eb 66/9-2) and the Austrian Science Fund 
(FWF grant Y328, START Prize) is gratefully acknowledged. 

References 

Abramowitz, M. and I. Stegun (Eds.) (1968). Handbook of Mathematical 
Functions (5th ed.). Dover. 

Applebaum, D. (2004). Levy Processes and Stochastic Calculus. Cam- 
bridge University Press. 

Barndorff-Nielsen, O. E. (1977). Exponentially decreasing distributions 
for the logarithm of particle size. Proc. R. Soc. Lond. A 353, 401-419. 

Barndorff-Nielsen, O. E. (1997). Normal inverse Gaussian distributions 
and stochastic volatility modelling. Scand. J. Statist. 24, 1-13. 

Barndorff-Nielsen, O. E., T. Mikosch, and S. Resnick (Eds.) (2001). Levy 
Processes: Theory and Applications. Birkhauser. 

Barndorff-Nielsen, O. E. and K. Prause (2001). Apparent scaling. Finance 
Stoch. 5, 103 - 113. 

Belomestny, D. and M. Reifi (2005). Optimal calibration for exponential 
Levy models. WIAS Preprint No. 1017. 

Bertoin, J. (1996). Levy processes. Cambridge University Press. 

Black, F. and M. Scholes (1973). The pricing of options and corporate 
liabilities. J. Polit. Econ. 81, 637-654. 

Borovkov, K. and A. Novikov (2002). On a new approach to calculating 
expectations for option pricing. J. Appl. Probab. 39, 889-895. 

Carr, P., H. Geman, D. B. Madan, and M. Yor (2002). The fine structure 
of asset returns: an empirical investigation. J. Business 75, 305-332. 

Carr, P. and D. B. Madan (1999). Option valuation using the fast Fourier 
transform. J. Comput. Finance 2(4), 61-73. 

Cont, R. (2001). Empirical properties of asset returns: stylized facts and 
statistical issues. Quant. Finance 1, 223-236. 

Cont, R. and P. Tankov (2003). Financial Modelling with Jump Processes. 
Chapman and Hall/CRC Press. 

Cont, R. and P. Tankov (2004). Nonparametric calibration of jump- 
diffusion option pricing models. J. Comput. Finance 7(3), 1-49. 

Cont, R. and P. Tankov (2006). Retrieving Levy processes from option 
prices: regularization of an ill-posed inverse problem. SIAM J. Control 
Optim. 45, 1-25. 

Corcuera, J. M., D. Nualart, and W. Schoutens (2005a). Completion of a 
Levy market by power-jump assets. Finance Stoch. 9, 109-127. 

Corcuera, J. M., D. Nualart, and W. Schoutens (2005b). Moment 
derivatives and Levy-type market completion. In A. Kyprianou, 
W. Schoutens, and P. Wilmott (Eds.), Exotic Option Pricing and Ad- 
vanced Levy Models, pp. 169-193. Wiley. 



48 ANTONIS PAPAPANTOLEON 

Devroye, L. (1986). Non-Uniform Random Variate Generation. Springer. 

Eberlein, E. (2001). Application of generalized hyperbolic Levy mo- 
tions to finance. In O. E. Barndorff-Nielsen, T. Mikosch, and S. I. 
Resnick (Eds.), Levy Processes: Theory and Applications, pp. 319-336. 
Birkhauser. 

Eberlein, E. (2007). Jump-type Levy processes. In T. G. Andersen, R. A. 
Davis, J. -P. Kreifi, and T. Mikosch (Eds.), Handbook of Financial Time 
Series. Springer, (forthcoming). 

Eberlein, E., K. Glau, and A. Papapantoleon (2008). Analysis of valuation 
formulae and applications to exotic options in Levy models. Preprint, 
TU Vienna (arXiv/0809.3405). 

Eberlein, E. and J. Jacod (1997). On the range of options prices. Finance 
Stoch. 1, 131-140. 

Eberlein, E. and U. Keller (1995). Hyperbolic distributions in finance. 
Bernoulli 1, 281-299. 

Eberlein, E. and F. Ozkan (2003). Time consistency of Levy models. 
Quant. Finance 3, 40-50. 

Eberlein, E. and A. Papapantoleon (2005). Symmetries and pricing of 
exotic options in Levy models. In A. Kyprianou, W. Schoutens, and 
P. Wilmott (Eds.), Exotic Option Pricing and Advanced Levy Models, 
pp. 99-128. Wiley. 

Eberlein, E., A. Papapantoleon, and A. N. Shiryaev (2008). On the duality 
principle in option pricing: semimartingale setting. Finance Stoch. 12, 
265-292. 

Eberlein, E. and K. Prause (2002). The generalized hyperbolic model: 
financial derivatives and risk measures. In H. Geman, D. Madan, 
S. Pliska, and T. Vorst (Eds.), Mathematical Finance - Bachelier Con- 
gress 2000, pp. 245-267. Springer. 

Eberlein, E. and E. A. v. Hammerstein (2004). Generalized hyperbolic 
and inverse Gaussian distributions: limiting cases and approximation 
of processes. In R. Dalang, M. Dozzi, and F. Russo (Eds.), Seminar 
on Stochastic Analysis, Random Fields and Applications IV, Progress 
in Probability 58, pp. 221-264. Birkhauser. 

Esche, F. and M. Schweizer (2005). Minimal entropy preserves the Levy 
property: how and why. Stochastic Process. Appl. 115, 299-327. 

Goll, T. and J. Kallsen (2000). Optimal portfolios for logarithmic utility. 
Stochastic Process. Appl. 89, 31-48. 

Goll, T. and L. Riischendorf (2001). Minimax and minimal distance mar- 
tingale measures and their relationship to portfolio optimization. Fi- 
nance Stoch. 5, 557-581. 

Hakala, J. and U. Wystup (2002). Heston's stochastic volatility model ap- 
plied to foreign exchange options. In J. Hakala and U. Wystup (Eds.), 
Foreign Exchange Risk, pp. 267-282. Risk Publications. 

Jacod, J. and A. N. Shiryaev (2003). Limit Theorems for Stochastic Pro- 
cesses (2nd ed.). Springer. 

Kallsen, J. (2000). Optimal portfolios for exponential Levy processes. 
Math. Meth. Oper. Res. 51, 357-374. 



INTRODUCTION TO LEVY PROCESSES 49 

Kallsen, J. and A. N. Shiryaev (2002). The cumulant process and Esscher's 
change of measure. Finance Stock. 6, 397-428. 

Kou, S. G. (2002). A jump diffusion model for option pricing. Manag. 
Sci. 48, 1086-1101. 

Kiichler, U. and S. Tappe (2008). Bilateral gamma distributions and pro- 
cesses in financial mathematics. Stochastic Process. Appl. 118, 261- 
283. 

Kyprianou, A. E. (2006). Introductory Lectures on Fluctuations of Levy 
Processes with Applications. Springer. 

Kyprianou, A. E., W. Schoutens, and P. Wilmott (Eds.) (2005). Exotic 
Option Pricing and Advanced Levy Models. Wiley. 

Madan, D. B. and E. Seneta (1990). The variance gamma (VG) model for 
share market returns. J. Business 63, 511-524. 

Matache, A.-M., C. Schwab, and T. P. Wihler (2005). Fast numerical 
solution of parabolic integro-differential equations with applications 
in finance. SIAM J. Sci. Comput. 27, 369-393. 

Matache, A.-M., T. v. Petersdorff, and C. Schwab (2004). Fast determin- 
istic pricing of options on Levy driven assets. M2AN Math. Model. 
Numer. Anal. 38, 37-72. 

Merton, R. C. (1973). Theory of rational option pricing. Bell J. Econ. 
Manag. Sci. 4, 141-183. 

Merton, R. C. (1976). Option pricing with discontinuous returns. Bell J. 
Financ. Econ. 3, 145-166. 

Papapantoleon, A. (2007). Applications of semimartingales and Levy pro- 
cesses in finance: duality and valuation. Ph. D. thesis, University of 
Freiburg. 

Prabhu, N. U. (1998). Stochastic Storage Processes (2nd ed.). Springer. 

Protter, P. (2004). Stochastic Integration and Differential Equations (3rd 
ed.). Springer. 

Raible, S. (2000). Levy processes in finance: theory, numerics, and empir- 
ical facts. Ph. D. thesis, University of Freiburg. 

Samorodnitsky, G. and M. Taqqu (1994). Stable non-Gaussian Random 
Processes. Chapman and Hall. 

Samuelson, P. A. (1965). Rational theory of warrant pricing. Indust. 
Manag. Rev. 6, 13-31. 

Sato, K. (1999). Levy Processes and Infinitely Divisible Distributions. 
Cambridge University Press. 

Schoutens, W. (2002). The Meixner process: theory and applications in 
finance. In O. E. Barndorff-Nielsen (Ed.), Mini-proceedings of the 2nd 
MaPhySto Conference on Levy Processes, pp. 237-241. 

Schoutens, W. (2003). Levy Processes in Finance: Pricing Financial 
Derivatives. Wiley. 

Schoutens, W. and J. L. Teugels (1998). Levy processes, polynomials and 
martingales. Coram. Statist. Stochastic Models 14, 335-349. 

Selivanov, A. V. (2005). On the martingale measures in exponential Levy 
models. Theory Probab. Appl. 49, 261-274. 

Shiryaev, A. N. (1999). Essentials of Stochastic Finance: Facts, Models, 
Theory. World Scientific. 



50 ANTONIS PAPAPANTOLEON 

Webber, N. (2005). Simulation methods with Levy processes. In A. Kypri- 
anou, W. Schoutens, and P. Wilmott (Eds.), Exotic Option Pricing and 
Advanced Levy Models, pp. 29-49. Wiley. 

Financial and Actuarial Mathematics, Vienna University of Technology, 
Wiedner Hauptstrasse 8/105, 1040 Vienna, Austria 
E-mail address: papapanOf am . tuwien . ac . at 
URL: http : //www . f am . tuwien . ac . at/~papapan 



